91/26 

February 1, 2008 



Non-linear WKB Analysis of the String Equation 



F. Fucito 

Dipartimento di Fisica, Universita di Roma II Tor Vergata and INFN, 
sezione di Roma Tor Vergata, Via Carnevale, 001 73 Roma, Italy 

A. Gamba 

Dottorato in Matematica, Universita di 
Milano, via Saldini 50, 20133 Milano, Italy 

M. Martellini^ 

INFN, sezione di Roma I, piazzale Aldo Moro 5, Roma, Italy 

O. Ragnisco* 

Dipartimento di Fisica, Universita di Roma I La Sapienza and 
INFN, sezione di Roma I, piazzale Aldo Moro 5, Roma, Italy 

ABSTRACT 

We apply non-linear WKB analysis to the study of the string equation. Even 
though the solutions obtained with this method are not exact, they approximate 
extremely well the true solutions, as we explicitly show using numerical simula- 
tions. "Physical" solutions are seen to be separatrices corresponding to degenerate 
Riemann surfaces. We obtain an analytic approximation in excellent agreement 
with the numerical solution found by Parisi et al. for the k = 3 case. 
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1. Introduction 

During the last two years, the partition function of two-dimensional gravity has 
been set in correspondence with the r-function of the KdV hierarchy, subjected 
to the constraint of the so-called string equation' 11 . In the case of pure gravity, 
described by a matrix model with criticality index k = 2, the string equation 
reduces to the well known Painleve type I equation for the specific heat of the 
theory. The string equation can be seen as a perturbation of a stationary KdV 

[2] 

equation, and thus solved in a semiclassical approximation around the stationary 

KdV solution. Novikov and Krichever' 31 have conjectured that such an approach 
lead to exact solutions already at 0-th order in perturbation theory. We have 
resorted to accurate numerical simulations in order to check this conjecture. On 
the other side, we show that "physical" solutions to the k-th multicritical model 
are obtained as separatrices corresponding to degenerate Riemann surfaces. We 
explicitly compare the k — 3 solution with the one found numerically by Parisi 
et al\ 4] . Purely asymptotic analysis seems to give in this case the same degree of 
precision of the semiclassical approximation in the non-degenerate case. 

1.1. The String Equation 
The string equation 

N 

-x + Y,U + i)tj+iRM*)] = o ( L1 ) 

3=0 

(where Rj are the Gel'fand-Dikii differential polynomials) for finite N gives a con- 
straint which is compatible with the first N flows of KdV. For N — > oo this gives 
a constraint compatible with all of the KdV flows. As a matter of fact, write the 
?-th KdV flow as Jp- — Kj[u(x)] = JfcRj[u>(x)], differentiate (1.1) with respect to 
x and compute 

8 N 

S 3=0 
N 

N 

sK s . 1 + K' s \Y^{j + l)t 3+ iK j } = 

sK a - 1 + K' a [-l]=0 
where K'A(f)](u) = j^Kj(u + e0)| e =o and we used known properties of the KdV 



1 



flows, in particular (i) commutativity, i.e. K'-[K S ] — K' s [Kj] = [Kj,K s ] = 0; (ii) 
K' g [l] = [K s ,l] = sK s -\, expressing the fact that 1 = r_i is the first master 
symmetry W of the KdV hierarchy. We note by the way that (1.2) can be written 

in a more satisfactory invariant form' 61 and easily generalized to arbitrary master- 
symmetries of the KdV equation. 

1.2. The Whitham method 

The Whitham method is the application to non-linear equations of the semi- 
classical approximation known in the realm of linear equations as WKB method. 
Let us consider e.g. the wave equation 

<Pxx - -tfPu = (1.3) 



and the plane wave solutions (p = Ae^ x u ^ . If we are studying propagation of 
light we have a "small" scale, i.e. the characteristic length of oscillation, which is 
about e = 10 -6 times smaller than the natural unit length used by the observer: 
so we can consider that solutions be locally given by plane waves, while on the 
"observational" scale the parameters A, k are "slowly" varying: 

,S(ex, et) 
i 

ip(x,t) = A(ex,et) e e (1.4) 

where |f = k(X,T), § = u(X,T), X = ex, T = et, £ = e&, | = e^. At 
0-th order, we find the eikonal equation of geometrical optics: (^f-) 2 = Jr(f^) 2 - 
The evolution of k = gives the paths of the "rays"; the eikonal equation is 
equivalent to 




1.5) 



In other words, we started with a class of exact solutions ip(x, t; A, k, u) = Ae^ x ^1 
of (1.3) and passed from the precise description of the oscillation process to the 
approximate description of the "slow" variation of parameters k,u given by (1.5); 
as a matter of fact we averaged over the rapid variation of the function (p and chose 
to observe only "secular" variations. The same method could have been readily 
applied to non-linear equations if we had 
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i) a family of exact solutions depending on an adequate number of parameters 
Ei, ... , E 2g +i; 

ii) a way of "averaging out" fast oscillations, in order to obtain an analog of the 
second equation in (1.5), which we will call "Whitham equation". 

This is the case for the KdV equation and more generally for equations that 
can be written in the form* 

___ + [LM ] = 0, (1.6) 

that is, as the compatibility condition for the existence of a solution ip of the linear 
system 



;i-7) 



where L = YJj=Q u j( x , V, A = Ylk=Q v k{ x ,y^)-^ are differential operators 

with scalar or matrix coefficients: as is known, these equations admit large sets of 
exact solutions (the so-called g-zone solutions) expressed in term of the Riemann 
6 function. 

In particular, the KdV equation 

4ut = u xxx — 6uu x (1.8) 
can be written as Lt = [L, A] with 
L = -d 2 + u(x,t) 

* ^3,^3,, d (1-9) 
A = d — —u(x, t)o — -u x (x, t), d= — 

and admits the "cnoidal wave" solution 

u(x,t) = p(x - vt;g 2 ,g3) (1.10) 

where p is the Weierstrass elliptic function and Q2,gs are arbitrary constants. 
Roughly speaking, the Whitham method consists in slowly varying the constants 



* [L, A] = LA — AL is the usual commutator of differential operators, while ^ and are 
supposed to act on the coefficients Uj , v/, ■ 
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g2,g3, that is in finding the correct dependence g<i = 92(X,T), g$ = gs(X,T), 
in order to approximate either (i) new solutions of (1.8) (corresponding e.g. to 
non-periodic initial data!) or (ii) solutions to the perturbed equation 



4u t = u xxx - 6uu x + eK (x) . (1-11) 

The cnoidal wave solution is a so-called 1-zone (g — 1) solution; general g-zone 
solutions are given by 



d 2 

u(x) = -2^-2 log 9(Ux\B) + C (1.12) 



where the exact form of U = V(Ej), B = B(Ej), C = C(Ej) is given in the 
Appendix. 



2. The Whitham method 

Let's start with perturbation theory for a non-linear equation (e.g. the sta- 
tionary KdV equation): 

u xxx - 12 uu x = eK (2.1) 
and look for solutions in the form 

u = uq + eu\ + e 2 «2 H (2.2) 

with 



u k = u k {t\X) = u k {^p-\E j {X)), X = ex, k = 0,1,2,... (2.3) 



where ^Si(X), . . . , ^^(X) are rapidly oscillating functions that will be determined 
in the following, and u k depends on some parameters Ej = Ej(X) which on their 
turn are slowly varying with x, and will be determined in the sequel. Substitute 
(2.2) and (2.3) in (2.1). Note that the form of the functions u k implies ^ = 
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^ • J^r + £-§x- Now take the various order in e: e.g. at order 0(1) we get 



(— --) 3 Mt|x)-i2Mt|x)( 



<9S 9 



) M0 (t|X) = 



(2.4) 



Let £ = -fix • = ~52j=i wt~5F~i an d S° on writing equations at all orders in the 
compact form: 



0(1) 
0(e) 
0(e fc ) 



£ 3 «o — 12 Mo £wo = 
£ 3 «i - 12 £(« «i) =Fi + K 
C\ k - 12 C(u u k ) =F k , k — 2,3,4, 



(2.5) 



Note that only the first equation is non-linear, and that equations for the order 
0(e fc ), k = 1,2,3,..., differ only in the non-homogenous term F k . This term 
would not be present in naive perturbation theory, and comes from differentiation 
with respect to slow variables: for instance uq x = (^ • + e^)«o produces 
the term of order e giving contribution to F\, together with derivatives like 
-§^S^ coming from uq xx , etc. (the explicit form of the term F\ is given in the 
Appendix). The idea is that at order 0(e) the term F\ should compensate for 
the perturbation term K in order that the "correction" u\ be bounded, and this 
gives equations for the correct dependence of the parameters Ej = Ej(X) on the 
slow variables; otherwise, we can think of averaging the 0(e) equation in (2.5) 
over the fast variables t\, . . . , t g , thus remaining with the only variable X. A third 
point of view is that the 0(e) equation has the form Lu\ — F\ + K, where L is a 
linear operator: this means that periodic solutions exist iff F\ + K is orthogonal to 
Ker(L^). The three points of view are all equivalent and give the same Whitham 
equations for Ej = Ej(X). 

Equation (2.4) has been written in explicit form in order to make clear an im- 
portant point: at any order the variable X appears as a parameter, while equations 
are only in the differential variables ti, . . . , t g . Moreover, putting 

BS 

— = XJ(E J (X)), t = Ux (2.6) 
we see that (2.4) becomes equivalent to the unperturbed uq xxx — I2u u 0x = 0, so 



* This explains why the variables x and X are usually treated as independent variables in the 
two-scale method, and makes rigorous the usual argument of "freezing" the slow variable. 
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we immediately get the form of uq^: 

«o(t) = -2(U- ^) 2 log^(t|^-) +C(E j ) (2.7) 

that is, in the genus g — 1 case, uo(*i) = fp(2to> • t\; g2,93)- 

Functions of the form (2.7) are periodic in t\, . . . ,t g ; moreover, we have seen 
that the variables t and X must be regarded as independent: so we can average 
both sides of the 0(e) equation in (2.5) over t\, . . . , t g , thus being left with the only 
variable X. The equations for E\(X), . . . , E2 g +i(X) thus obtained in the case of 
the stationary KdV equation have the form 

8 , E+r(E 1 (X),E 2 (X),E 3 (X)) 

— J(E-E 1 {X)){E-E 2 {X)){E-E z {X))dE = -6 - dE 
8X W ^/(E-E 1 (X))(E-E 2 (X))(E-E 3 (X)) 

(2.8) 

where r(E u E 2 , E 3 ) = /f 2 , EdE / ff 2 , dE and 

E is a dummy variable. 

In general, the Whitham equations for KdV are differential conditions on the 
functions p = p(E, Ej(X,T)), Q = Q(E, Ej(X,T)) (quasi-momentum and quasi- 

energy) appearing in the E — > oo leading term ip ~ eP x ^ of the solutions to 
the associated linear system (1.7). For KdV they take the form 

dfl dp _ (^Kip) dp 

8X 8T (flij,) dE { ' } 

Equation (2.8) can equivalently be written as ^ = —6^. A complete derivation 
and explanation of these formulas is given in the Appendix. 



f Actually, this is the complex form of uq. A rigorous treatment needs the use of real variables 
ti, . . . , t 2g - In the 3 = 1 case we would get uo{t\,t2) = p{2uj ■ t\ + 2a/ • t 2 ; 32, 33)- 
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3. The Painleve type I equation 



The Painleve equation 

u" = 6(u 2 -x) (3.1) 

(which is obtained from (1.1) when we leave only ts ^ 0) can be differentiated once 
with respect to x and seen as a perturbation of the (integrable) stationary KdV 
equation: 

u'" = 12m/ + eK (3.2) 

with e = 1, K — —6. An asymptotic relation between the two equations is obtained 
through the change of variables 



It comes out(2.9) t 181 that solutions of (3.1) are asymptotic to functions of the 
form u(x) = p(|x 5 / 4 ; 12, (73); nothing is said about the parameter g%. We will now 
resort to a finer analysis of the problem. 

We are here concerned(2.9) | [9] with the case of the torus of equation (2w) 2 = 
4E 3 — giE — gs = 4(E — E\)(E — E2)(E — E3); we distinguish two cases; for 
A = g\ — 27 g\ > we have real roots E\,E2, £3 and two periods 2oo and 2uo' which 
are respectively real and pure imaginary (2. 9) | [101 : 



E2 Ei 

f dE 2K(m) , f dE 2iK(l-m) 

2u = 2 - = ' . 2J = 2 \ - = 

J ^4:E 3 -g 2 E-g 3 JW-E 3 J y/±E s -g 2 E-g 3 VE1-E3 

(3.4) 

where 



K{m) = / — — — , E(m) — / y 1 — msin (f)d(j), m 



l-m sin 2 6 J E\ - E s 

o v o 

(3.5) 

are standard elliptic integrals of the 1 st and 2 nd kind, respectively, and m is called 
the Jacobi modulus. We remind that the solutions of (3.2) with e = can be 



expressed in term of the Weierstrass elliptic function (see Ref. 10) 

P^92,93) = -^lo g e l (^-)-^. (3.6) 

The Whitham equations for (3.2) take the form 

^dE = -6dp (3.7) 

or, explicitly, 

Oy/uP-nE-gsdE = dX m dE = -6 - E + r(X) dE 
dX - ' _ g2 E - g 3 y/4E 3 - g 2 E - g 3 

(3.8) 

giving 

|| = 12, |f = 12 r (X). (3.9) 

The first equation is readily integrated, giving g 2 = 12 • X + const; moreover, 
the second too is integrable by quadratures, because §^dE already has the same 
behaviour as — 6dp for E — > oo, so we need only to impose the normalization 
condition (compare with (5.4)) 

Im^jwdE^j = Re(j J - 12XE - g 3 (X) dE^j = const = h (3.10) 

in order to get Im § ^dE = ^Im § wdE = and ^wdE = -6dp. But (3.10) 
can be solved explicitly for g%, yielding the correct dependence g% = g%(X). 

In order to do this it is convenient to introduce new parameters (A, m) in place 
of (92,93)- 

(Ei — E3) 2 E2 — £3 

A = , m = — — , 

9 E\ — E$ 

g 2 = 12A(1 -m + m 2 ), g 3 = 4A 3/2 (2 - 3m - 3m 2 + 2m 3 ), 

Ei = v / A(2-m), E 2 = v / A(-l + 2m), E 3 = V\(-l - m), (3.11) 

Ei 

r A 5 / 4 
J y/AE* - g 2 E -g 3 dE = i—^ m )- 

E2 



giving (for A > 0) the solution 



*<•»> = I |g I • (3-12) 



In Fig. 1 we show the form of the resulting function x = x{m), both for the cases 
A > and A < (for graphic convenience we plotted x versus ^ instead of m). 
The pole in the A < region gives rise to two distinct curves in the space of 
parameters (<72,<73), which we plotted in Fig. 2. Knowing the properties(2.9) f [111 
of the Weierstrass elliptic function p we see that for x = || > the Whitham 
method gives us either small oscillations (curve in the A > region) or a sequence 
of poles (lower part of the curve, lying in the A < region). If we start with 
oscillatory behaviour at +oo and go in the direction of decreasing x, we ultimately 
reach A = where a transition to polar behaviour occurs. Note that for x < we 
can have only polar behaviour. The two parts of the curve correspond to solutions 
having or not having poles at +oo. The "physical" solutions will be seen to be 
separatrices lying between these two kinds of solutions. 

The function $(m) is computed reducing the integral in (3.10) to standard 
elliptic integrals(2.9) T™. We get /| 2 g = *$L, g = K(l-m)y/3\V\ 

and finally 



^2 _ , . 

/2 A / 
^4^3 - g 2 E - g 3 = -——(18(2 - 3m + m 2 )K(m) - 36(1 -m + m 2 )E(m)), 
5 v3 

Es 
E 1 



/2 iA 5 / 4 
- g 2 E - g 3 = -—^(-18(m + m 2 )K(l - m) + 36(1 -m + m 2 )E{l - m)), 
5 V3 

E'z 



$( m ) = -18(m + m 2 )K{l - m) + 36(1 - m + m 2 )E{l - m). 

(3.13) 

This procedure gives us also the "actions" S\,S2'- 



Ei E 2 
Si = ~J V / ^-g 2 E-g 3 dE, S 2 = ~ J y/AE^ - g 2 E - g 3 dE, (3.14) 

E'2 E3 
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because from jj^dE = —6dp it follows 



<)Sl Ui = 2 f dp, ^ ( ■> 2 I dp. (3.15) 



dX J " dX 

E 2 E A 

Take S = S\ + tS 2 as prescribed by (2.6), (5.15) and (5.16), where r = i K ^ r ^ > 

is the 1-dimensional analog of the period matrix B, and find from (3.13) and the 
Legendre relation EK' + E'K - KK' = § (see Ref. 10): 

^ B/4 .12(l-m + m») = 4 x 
5 K(m) 5 2w y J 

The Whitham solution has the form 

,4 ._, , _ 4 



M 



V) = + cu; 12a;, 4/(:r):r 3/2 ) = -y/xp(-ix 5 / 4 + uj- 12, 4J(a;)); (3.17) 



(We have used here the homogeneity property p(tx; p-, ||) = -pp(x; g 2 , #3), t = 

rr -1 / 4 ; the upshift is needed in order to get non-singular solutions). We thus recover 
the "Boutroux" asymptotic form (3.3), but with a more precise 

w / u 2-3m-3m 2 + 2m 3 

I(x) = I (m(x)) = — /2 ■ (3.18) 

(1 — m + m 1 ) 6 ! 1 

Asymptotic analysis(2.9) | [131 shows that the amplitude of the oscillatory solu- 
tions of (3.1) decreases at +00 as -4/s- This feature, not recovered by the simple 
Boutroux-type solution, is obtained from the fine tuning realized by the term I(x). 
Just use the expansion 

K{m) = |(1 + ^ + ^m 2 + ■ ■ ■), E(m) = 1(1-^ - + . . •), (3.19) 

for m — > (see Ref. 10). From (3.11), in the limit to — > 1, we get x = || ~ A(m), 
$(to) ~ ±f^(l - to) 2 , A(to) ~ (%^) 4 / 5 (l - to)~ 8 / 5 , -E u -E 2 ~ Jx,E x -E 2 ~ 
\J ~7^ x ~ 1 ^ (Notice that the function p(ix + a;) oscillates between the extremal 
values —Ei and —E 2 , see Ref. 11). 
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If we want to examine the case x < 0, corresponding to A < 0, formulas (3.11) 
are no longer convenient, and we resort to the following real parameterization 



m* = — = \ + ia, A* = m 2 A = -p 2 , iA 5/4 $(m) = -iA* /4 $(™*) (3.20) 
m 2 

The solution corresponding to (3.12) in the complex case becomes 



4/5 



^ = " 1 ' (3 ' 21) 

A numerical computation shows that the function $(m*) = $(| = $(c) has 
a zero for a = <7 = -0.231026398427 .... 

In the case A < real and pure imaginary combinations of periods are given 

Kv / , — ^M) , / - iK(l-rh) h rx2 _ q P 2 £2 _ 1 3-4a 2 i ~ _ 1 , 



/ff ' 2 v/H ' 1 4 '1+4^1'"" 2 1 yTT4^2 

(see Ref. 10). With this parameterization we recover the numerical result of 
Ref. 13, asserting that the distance of poles goes asymptotically as -t/i, with c = 



7.276726 . . .; as a matter of fact we find c = = \K(\ - / a ° J(|^) 1/4 = 

2.970711275212..., which exactly coincides with the result of Ref. 13 after the 
rescaling c = \/6c (due to our factor 6 in (3.1)). An analogous reasoning for x > 
gives the period of the oscillatory solutions going as ^74, with c—^j- — ^. 

In Fig. 3 we show the approximate solution (with h — 1), together with an 
exact numerical solution obtained with the Runge-Kutta method (the Painleve 
equation is satisfied with an error of 10~ 14 ). Fig. 4 is a magnification of the region 
around the zero where the approximation seems to be less effective. We see that 
the solution we are considering seems to be out of phase with respect to the exact 
numerical one. This is no surprise because we really have neglected a phase: the 
Whitham method made variating the "constants" #2 and #3, but in p{x + c; gi, #3) 
we have also a third integration constant c, that we assumed to be zero. We do 
believe that an equation for this phase can be deduced(2.9) | 1141 , and its behaviour 
will be discussed in a forthcoming paper. Here we have just made a fit of such a 
phase on the side of the positive x, getting a correction 8(f>(x) ~ j— ^^pgg ■ We 
have also done the same on the negative x side but we got a really tiny correction 
that we chose to neglect. The fit of the "experimental" data is shown in Fig. 5, 
while Fig. 6 shows the effect of putting in the correction by hand (for x ~ we let 
6<f>(x) die smoothly). 



11 



4. Degenerate solutions 



The solutions u(x) of Painleve-like equations represent specific heats of the ran- 
dom matrix models, which in the planar limit must satisfy the boundary condition 
u(x) ~ x 1 ^ (x — > oo) for scaling arguments of the partition function: Z ~ x~ 1+2 , 
where 7 is the string susceptibility. Thus u(x) = 0^ ~ -§^ x ~ 1+1 ~ x ~ 7 > an d the 
result comes from the fact that 7 = -| in the proximity of the critical point. The 
situation is common to other non-linear physical models in the critical regime. It 
is well-known that the problems of mathematical physics must be complemented 
by boundary conditions, and that the boundary conditions contain in some sense 
the physics of the problem. In our particular case we come to the request that the 
solutions to the string equation for the k-th multicritical model (which is obtained 
from (1.1) putting all tj = 0, except tfc+i) must satisfy the physical constraint 
u(x) ~ x 1 ^ for x — > +00. However, the Whitham method gave us either small 
oscillations modulating over — y/x (see Fig. 3), or solutions with poles. The only 
possibility to get non-periodic solutions is to considerate degenerate Riemann sur- 
faces, where the length of the bands is sent to 0. These solutions are degenerate 
cases of the periodic solutions and are themselves unstable separatrices, lying be- 
tween the two sets of solutions with poles and without poles for x — > +00. 

All we have to do is to compute the spectral curve corresponding to the given 
k-th stationary KdV equation and imposing the coincidence of the pairs of branch 
points; g — k — 1 conditions are found by requesting that dw = dE 9+ 2 + (9(1); one 
more condition comes from fixing the periods of the solutions at x = ±00; the last 
g conditions come from the request that the branch points coalesce in pairs. 

For the k = 2 case we require dw ~ dE 3 / 2 + 0(1), giving 

Ei + E2 + E% = ,^ 
E1E2 + E2E3 + E3E1 = c. 

The asymptotic condition fixes c = 3, and we ask for E2 = -E3. This gives E2 = 
E3 = — -j- = 1. The 6 function degenerate to a combination of hyperbolic functions 
and we finally get 

(Vx~(l — 5 — ), x>0, 

{y/=x~(p(l(-x)-4; -12, -I(± + ia )), x < 0. 

(In the x < case the surface does not degenerate, but I(x) — > J(| -Mo"o) — const, 
where 00 is the constant introduced in the previous section). The function u(x) is 
plotted in Fig. 7. In x = we get a cusp as we are trying to connect at finite x 
two asymptotic solutions: a smoother curve would probably require a phase-type 
correction as suggested in the previous section. 
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In the k = 3 case the condition dw = dE^I 2 + 0(1) gives 



fEti^ = o 

\ J2i<j EiEj = 

( Y.i<j<k E i E j E k = C 

Asymptotic conditions fix c = ±| for x — > ±oo; we put E<i = E% = s, E% = E4 = t, 

and find s = ±( 1+ ^ ), t = s, i?i = ±1 (rr — > ±00). Again the 9 function factorizes 
in products of trigonometric and hyperbolic functions (for details see(2.9) |' 151 
and(2.9) | 1161 ) giving 



Vq(x) 



2 ( — ) sinh(a;r+21n ^) sin(fex)+2 cosh(ax+2 In ^) cos{bx)— 2 

~ (sinh(aa; + 21nf)-f sin(6x))2 ' X < 0; 

1 26^ ^ sinh(6;r+21n ^) sin(aa;)+2 cosh(&E+2 In ^) cos(ax)— 2 
(Sinh(frr + 2 ln^)-^sin(aa;)) 2 



(4.3) 

where a = (30)2 cos |, 6 = (30)2 sin |, ^ = arctan-^g. (We checked that this so- 
lution satisfies the higher order stationary KdV equation). Fig. 8 shows that 
uq(x) = 3 y / a; • vq(^x&) approximate very well the form of the solution found nu- 
merically in Ref. 4, except that in the proximity of x — 0. (In the graph we have 
shifted z\ \— > z\ + \, Z2 i— > Z2 — \ for x < to match the phase of Ref 4. The 
resulting function still satisfies the KdV equation). 
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5. Appendix 

In sections 1, 2, 3 we review some facts about algebraic geometry and KdV 
equations, mainly for notational convenience. In section 4 we report the proof of 
Krichever's theorem, following Ref. 2. 
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5.1. Complex curves 
The algebraic equation 

w 2 = E 29+1 + aiE 2 ^ 1 + a 2 E 29 ~ 2 + ... + a 2g -iE + a 2g = P2 g +i (E) (5.1) 

defines a curve T in the complex plane of the variables (E,w). The curve is com- 
pactified at the oo and is known to be topologically equivalent to a compact surface 
with g holes. If the polynomial p2 g +i(E) has 2g + l real roots Ei, ... , E2 g +i, we can 
draw them on the complex plane and use solid lines for the segments (E 2 k-i, E 2 k), 
where the square root w = ±y/p2 g +i(E) = ±\J(E - E\)(E - E 2 ) ■ ■ ■ (E — -E^+i) 
takes real values. In spectral theory these are the forbidden zones of the spectrum. 

Coordinates on T are given by 

almost everywhere, 

in the neighborhood of Ej, (52) 
in the neighborhood of 00. 

Consider integrals of the form f ' Vtu = f — . = dE , k = 0, 1, 2, . . . Us- 

6 2y/(E-E 1 )...(E-E 2g +iy 

ing (5.2) it is easy to see that Qq, . . . ,Q g ^i are everywhere non-singular, while 
fl g , flg+i, . . . have poles at the infinity of order 2, 4, . . ., etc. 

To fix a basis of differentials we chose first a canonical(2.9) | 1171 basis of 
paths a\, . . . , a g and b\, . . . , b g on T and take u>i, . . . , uj g as linear combinations 
of Qq, . . . , fl g -\ satisfying the normalization condition 

E^k 

ujj = 2 J u k = S jk , j, k = 1, . . . , g (5.3) 

a k E 2 k-i 

We are left with the 6-periods, forming a g x g matrix Bj^ = B k j = § hk ooj = 

2 f^ +1 Wj. Differentials with poles of order 2j will be indicated by u/ 2j ' -1 ) for 
future commodity, and can be fixed by requiring that they go at the infinity as 
u; {2j-i) ^ dEi~2+0(l). the arbitrariness on the holomorphic tail can be eliminated 
by imposing the 2g real conditions 



/ CJ (2J- 1 ) = 0, Im (fj 2j -^ =0, k = l,...,g. (5.4) 



Im 

a k b k 

(Another standard choice of the normalization is to impose instead of (5.4) the g 
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complex conditions ^ = 0, k — 1, . . . , g). We will also use the notation 

dp = cjW ~ dy/E, dVL = ~ dE 3 / 2 , E -> oo; (5.5) 

these are the differentials of the quasi-momentum p(E) and quasi-energy £l(E), 

fundamental in the theory of the KdV equation (see 5.11); p(E) = f^dp and 

Q(E) = f^dtt are multivalued functions with periods Uj = § h . dp, Wj = § h dQ, 

j — 1, . . . , g. Note that p(E) and Q(E) are uniquely determined by the asymptotic 
behaviour and the normalization conditions. 

5.2. Functions on the surface r 

The Abel map P i— > A(P) = t (J P . . . , f P ujg) maps any point P on Y on 
the g-dimensional torus C 5 /{period lattice}. The Abel map can be inverted by 
means of the Fourier series 

9(x\B) = J2 e nin - B - n+27vin - z (5.6) 
neC 9 

defining (for positive definite ImB) the Riemann function, which has the period- 
icity properties 0(z + ej) = 0(z), 0(z + Bej) = e _,riB "- 2,ri ^0(z) (see Ref. 17). As 
a matter of fact, a theorem of Jacobi asserts that the function 



f(P;P h ...,P g )=9( u k -J2 ^k + K k ) (5.7) 

1=1 

OO J A oo 

(where K = (K k ) is a certain constant vector) has exactly g zeroes P\,...,P g . 
This theorem gives an analog of the development of a rational function in simple 
fractions. 

5.3. The KdV equation 

The KdV equation 4ut = u xxx — 6uu x admits the Lax representation Lt = 
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[L, A], with L, A, given by (1.9), and has exact solutions of the form 



d 2 

u{x,t) = -2 

where (see (5.5)) 



u(x,t) = -2-^2 \og6(Ux + Wt + z |P) + C(B) (5.8) 



Uj=<bdp, W j= (bdn, j = l,...,g. (5.9) 



bj bj 

This can be seen as follows(2.9) j [18] . The equation Lt = [L, A] is the compatibility 
condition for the system of linear equations 

Lip = Eip 

M, = * < 5 - 10 > 

at 

For E ~ oo we get L ~ d 2 , A ~ <9 3 , so the asymptotic form of the common 
eigenvectors ip will be 

ip{x,t;E) =eP^ x + Q ^ t -(f)(x,t;E), with p(E)^v^E, {1(E) ~ 22 3 / 2 as £ 

(5.11) 

The exact form of (p(x,t) is (compare with (5.7); see Ref. 9) 

, gag - e; = i (Pj) + ^ + ^ + g*) ,_ 19 . 

0(x, t; 22) = ~fT~p~\ ( 5 - 12 ) 

E^iJ* (P 'W**) 

for given T and Pi, . . . , P g on T. It is easy to verify that ip = eP x ^ • <p is a 
one-valued function of 22. The function ip is uniquely determined by the behaviour 
at infinity (ip ~ e s 7 :r+£; 7 *) and the position of the g poles P\, . . . , P g , as can be 
easily seen with the help of the Riemann-Roch theorem. For E = k 2 — > oo we get 

iKs, t; 22) = ce** + kh • (1 + *M + + 

/e Ar 

If ^ satisfies = Eip, Aip = then we can collect terms of the same order in ^ 
and £i,£2i • • • should satisfy some equation at any order. It is easily seen that the 
first of these equation gives u(x, t) = 2^(x,t), so for this choice of the potential 
u we get (5.10) verified at order O(^). But note that (L — E)ip, (A — ^)ip again 

have the behaviour e kx + ^ ^ and poles at Pi, ... , P g , so they must again have the 
form (5.13) with no 0(1) term: this means that they are identically zero; so (5.11) 
and (5.12) give an exact solution to (5.10). Developing (5.12) at 1 st order in ^ we 
get precisely formula (5.8) for u = 2J^£i(x,t). 



oo.i 
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5.4. Krichever's theorem 



We will here report the form found by Krichever for the Whitham equations 
of systems of KP type, following Ref. 2. It is convenient to consider directly 
the general case of equations of the form (1.6), as for instance the KP equation 
Suyy + -^{Aut — 6uu x + u xxx ) = 0. (1.6) is the compatibility condition for the 
existence of a solution ip of the linear system Lip = Aip = The common 
eigenvector tp will be given here using a particular real normalization, necessary 
for the successive averaging procedure: 

ip(x,y,t;P) = eP x + E y + nt + s ' t -(j)(Ux + Vy + Wt + t,P) (5.14) 

here the spectral curve T is no more hyperelliptic, and consequently the spectral 
parameter E becomes itself a multi- valued function E(P) of the point P on the 
surface; the functions p(P), E(P), £l(P) are normalized by requiring that they have 
pure imaginary periods along all cycles a\, . . . , a g , b\, . . . , b g ; U, V, W are the real 
2g-dimensional vectors of periods of the multi- valued functions p(P), E(P), Q(P): 

U — t (^ dp, . . . , <j> dp, — <j> dp, . . . , — <j> dp), etc.; (5.15) 

b\ b g ai a g 

t\, . . . ,t2g are auxiliary "times" needed for further procedure of averaging; s\, . . . , S2g| 
are the corresponding "momenta" , not needed in what follows; is a periodic func- 
tion with period 1 with respect to all of the 2g variables x\, . . . , x g , y\, . . . , y g (A 
is the Abel map): 



0( 



(A. p) _ rp 2niA(P)-y g(A(P) + x + By + z )6(z ) 

Vy)' ] '0(A(P)+zo)0(x + 5y + z o ) ^ 



We will also need the solutions to the adjoint system ip^L = —^gy~, ift A = 
where differential operators written on the left should be intended according to 
^(u-^j) = (— -^y(ip^u). (formal integration by parts). The left and right action 
differ only for a complete derivative: 

ft f$ 
{^L)iP = ^LiP + _(^t L (ity) + _(^t L (2)^) + . . . ( 517 ) 

as can readily be seen by repeated applications of the Leibnitz rule. Here = 
" ri d(d) r ^ (f° rma l derivation with respect to the symbol d: for instance, A^> = 
-3d 2 + \u). 
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Solutions to the adjoint system can be written in the form ij;*(x,y,t; P) = 
e -px-Ey-m-s-t . 0t(_Ux -Vy-Wt-t, P) 

The application of the Whitham method to equations of the KP type is allowed 
by the possibility of averaging identities by means of some "ergodic theorem" : the 
average of periodic functions <f>(ti, . . . , t2 g ; Ik) with period 1 with respect to all 
arguments is given by (0) = / (p(t)d 2 n, and for a generic vector U it coincides with 
the limit {4>) x = lim Xo -»+oo I-° 0(Ux)da;, because the line Ux winds densely 
on a 2g-dimensional torus. Note that the derivative of along any direction has 
zero average value: 



2 g i 

{^(Vx)) = U • <^> = £ U S J ^d tl . . . dt 2g = (5.18) 



Note also that the functions ^j- are again periodic, because in our real normaliza- 
tion the periods have the fixed value 1 (not depending on the 

Solutions of the KP equation have the form 

d 2 

u(x,y,t) = -2^-2 log 0(U x + Vy + Wt + z \B) + C(B) (5.19) 
Thus, we look for a semiclassical approximation in the form 

u (x,y,t) = -2 log 6( S(X ' e F ' T) \I(X, Y, T)) + C(X, Y, T), 

where || = U, ff- = V, f| = W, = e^ s ' s 0(|). The operators L , A obtained 
substituting u \— > uq are taken as first terms of the asymptotic series 

A = A + eAi + ■ ■ ■ , L = L + eLi + ■ ■ ■ , (5.20) 



Let us introduce the notation 



<9r <9r (91 ^ <9r 97, 

J J 



then the substitution u h- > mq implies ^ h- > J| + taking terms of order 0(e) 
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in (1.6) we get the linearized equation 



w ~ it + Al] + [Ll ' Ao] = K ~ F ' (5 - 22) 

where F is the term due to derivation with respect to slow variables: this term 
must be adjusted in order to compensate for K. We readily find for F the form 



F= dL_dA {l) dA_ A(1) dL. 
8T BY y dX dX' 



(5.23) 



Now use % = Lip, % = -^L, f = °& = —i^A, and (5.17): 
|(^)-^(^) = ^ t (^-^ + ^o,A 1 ] + [L 1 ,A ])^+^(...) (5.24) 



Thus, the average of the left hand side of (5.22) comes out to be zero (being the 
average of a total derivative) and we obtain the Whitham equations in the implicit 
form 

l^Kij)) = (^F4>) (5.25) 

Explicit computing of the Whitham term (ip^Fip) will give us the final form. Take 
respectively 

(i) a curve I = I(r) in the space of parameters (P = const, t = const); 

(ii) a curve P = P{r) moving the point P on the surface T (I = const, t = 
const); 

(in) a curve t = t(r) moving only the "times" U (I = const, P = const). 
Correspondingly, we get L(t),A(t), etc., and 

^(r) = A T ) X + E ( T ^ + W + s • * • (P(IJ(t)x + Y(r)y + W(r)t + t) 
^ = e -px-Ey-Qt-s-t . (f) ^_ lJx _ Yy _ Wt _ t) 

Now compute J^| r=0 ip^ip(r) in all the three cases and use a point to denote deriva- 
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tion with respect to r: 



n 9 



^ty( r ) = (p x + E y + (it) ^ + (Ux + Vy + Wt) ■ + i • V^Sp 



in. 



d_ 



l-\ 9 



^ty( r ) = (dp x + dE y + dVt t) $(f>, 

i 

^( r ) =S .t0t + t^.t; 

, at 



Using = A^, fgj - and (5.17), see that 



(5.27) 







^tyty(r)) = ^(A(r) - Aty(r) - AtyU(ity( T )) + |_(. . .) ( 5 .28) 
Deriving the left-hand side with respect to r and using (i) we get for instance 



d 2 



d_ 



d 



^(^V(r))=QVV + W-V t ^ 



dip 



dt 

d 

+ {(px + Ey + ttt) — (^ty) ^ 5 - 29 ^ 



+ (Ux + Vy + Wt)4(^)+ i 4(^)> 



Now fix x, y, t and average upon rf 25 t: the terms in braces are linear combinations 
of total derivatives with constant coefficients and thus vanish, giving 



k 8t 



(5.30) 



this passage contains the essence of the method of averaging. We can now go on 
deriving both sides of (5.28) with respect to r, applying (i), (ii), (Hi) and finding 
after averaging the following identities: 



(i) : fltyty) + w • (^> = (^|^> - p^Wty) - u • <0 f i (1) ^> 

(ii) : dfl(^ip) = -dp{^A W ip) 
(Hi) : = <^> 

(5.31) 

Note that ip^ip = f 0, that we posed = e (~P x ),4We(P x + " '), and that 
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(Hi) implies 

r)A r)A 8 A ft A 

<^t|_^) = (^(-—A + (Us + Vy + Wt) ■ °—)^) = <^t|_^) (5.32) 

We get analogous identities for E,\,L if we derive instead of (5.28) the identity 
!^ty(r)) = ^(L(r) - L)^(r) - JV^^t)) + J^(. . .) (5.33) 

Now letting r = Y, T (remember that Y, T are independent of y, t), we can rewrite 
(5.31) and the analogous identities for L as 

(5.34) 

The last term we need comes from the identity 



d_ 

dt 



(^L^4,(t))-^A^4,(t)) = ^[L^\A(r)-A)-A^(L(r)-L)]i,(r)+^-(. . )| 

(5.35) 

which, after putting r = X and averaging, gives 



(5.36) 



Summing up (5.34) and (5.36), and using the compatibility conditions 
<9U dW dU dW d\ dW 



dY dX ' dT OX ' dT dY 

we get 



(5.37) 



= («-^)(^> + (« - ^L<%)H % -§WU<%,) (5.38) 

using (ii) from (5.31) we can rewrite the Whitham equations in the final form 

,dVt dE.j .dp dQ, Jrn ,dE <9p. JO (^K^) J 

" + " ^ + (^ - = (5.39) 
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5.5. Whitham equations for the stationary and evolutive KdV 



The KdV is a particular case of the KP equation. The solutions of the KP 
equation correspond to generic (non-hyperelliptic, that is not of the form w 2 = 
p(E)) Riemann surfaces T: in this case the function E(P),P G T itself is no 
more one- valued and its differential dE has non-zero periods Vj = ^ § h dE. The 

solutions have the form (5.19). When V is hyperelliptic the differential dE is exact, 
V = and the dependence on y disappears, giving solutions to the KdV equation. 
For KdV the Whitham equations have the form 

f dp dQ _ {^Kjj} dp 

[ dT dX>~ (V^> dE 15 UJ 

Solutions to the stationary KdV equation, which is equivalent to the linear system 
Lip = Eip, Aip = w(E)ip, come out when dQ too is exact, and this is true for 
Q(E) = 2w(E)dE = -\f AE 3 — g2E — #3 dE] the corresponding Whitham equation 
gives 

dw (tp^Kip) 

ax dE= -WW dp (5 ' 41) 

where dp is normalized with Im^ dp = 0. 
5.6. Integrability of the Whitham equations for K=0 

For KdV the Whitham equations have the form (5.40). In this case, the param- 
eters Ij of the preceding section are simply the branch points E = t (E\, . . . , i?2 5 +i) 
of the spectral curve. For K = this comes out as 

dp 

df = dx (5 ' 42) 

Krichever has shown that (5.42) has solutions E\(X, T), . . . , E2 g +i(X, T) given 
implicitly by the conditions 

^(E(X,T))U =S . (X;T) +X+T^(E(X,T))U =£ . (X)T) = 0, 3 = l,...,2i/+l, 

P P (5.43) 

where dA is an arbitrary differential with possibly discontinuities and singularities 
not depending on X, T. Analogous solutions exist for the KP case (see Ref. 2). 
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In order to see it, consider that if the function S(X,T) = f**dA(X,T) + 
Xdp(X,T)+Tdtt(X,T) is such that |f = p, ff = SI, then(5.42) is automatically 
satisfied. Now, (5.43) is equivalent to dS\ E=E .t x ,T) — f° r J ' — 1? • • • > 2g + 1; 
the form -^dS = dp + (-^dA + X-^dp + T-J^dQ) has the same normalization 
as dp, the same singularity at the infinity and is holomorphic everywhere else, ex- 
cept, generally speaking, in the points Ej, because e.g. ^( E ~e 1 )(e E e 2 )(e-e 3 ) = 

~ \ {(E-Ei)(E-E2)(E-E 3 )] 3 / 2 + ' ' ' cease to be holomorphic in these points. However, 
thanks to dS\ E =E j = we have dS = y/E - EjdE + 0(E - Ej) 3 / 2 in a neigh- 
borhood of Ej, and ^dS = 2 J E _ E +■■■ = -jrxS^ comes out to be 
non-singular in Ej too: thus -^dS coincides with dp, and similarly J^pdS = dfl. 

5.7. Integrability of the stationary Whitham equation for K — 1 

The equation §^dE = dp can be integrated in the following way (see Ref. 
3): take k = y/E, dE = 2k dk, w = y/ ' E 2 9 +1 + ciE 2g -i + ■ ■ ■ + c 2g = k 2 ^ +1 + 
2 S+lT 2g+ ik 2 ^ 1 + 2 9 T lT 2 g-ik 2 3' :i +. ■ .+ 3 T 3 k+^+0(^) and require Im/ 0i wdE =| 
hi = const, \m§ b ^w dE = h[ = const, i = l,...,g. The first condition gives 
ci, . . . , c g as algebraic functions of T 2g+ \, . . . , T3; the second condition fixes also 
Cg+i, . . . , c 2g as transcendental functions of T 2g+ i, . . . , T3, X. Moreover, the first 
condition gives -^wdE ~ d\f~E ~ dp, while the second assures that both sides 
have the same normalization: this means that -^wdE = dp, as desired. (Intro- 
ducing the functions Qj ~ A/ 7 + (9(1), E — > 00, it is easy to see that Jriy dE = dQj, 

thus giving also ^p 1 = ^f). 
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